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DISPERSED FOURIER TRANSFORM SPECTROMETER 

The present application claims the benefit of the priority filing date of provisional 
patent application no. 60/435,730, filed December 20, 2002, incorporated herein by 
10 reference. 

Field of the Invention 

This invention relates to a method and device for Fourier transform spectrometry. 
1 5 More particularly, the invention relates to a spectrally dispersed Fourier Transform 
Spectrometer. 

Background of the Invention 

20 The technique of spectrometry is used widely to determine spectra either occurring 

in nature or in laboratory settings. Recent advances have provided significant 
improvements in spectrometric applications such as astronomical spectrometry. Prior to 
1980, measurements of Doppler velocity shifts providing velocity precisions on the order 
of 1 km s" 1 were seldom possible. Now, using high precision absorption-cell 

25 spectrometers, measurement of velocities with precisions as small as 1 m s" 1 are attainted. 
These data make possible the detection of planetary companions to stars on the order of 
0.1 6Mj < Msin(i) < 15Mj (Mj is the mass of Jupiter, and i is the inclination angle of the 
orbit of the planetary companion). It would be desirable, however, to obtain a greater 
sensitivity capability to detect smaller companions or to provide greater sensitivity in other 

30 types of non-astronomical applications. 

Another available technique is Fourier Transform Spectrometry (FTS). An FTS 
spectrometer is an autocorrelation, or time-domain, interferometer. The theoretical basis 
was laid at the end of the 19th century (Michelson, A.A. 1891, Phil. Mag., 31, 256, 
Michelson, A.A. 1892, Phil. Mag., 34, 280.), but FTSs did not achieve widespread use 

35 until approximately 75 years later (Brault, J.W. 1985, in High Resolution in Astronomy, 
15th Advanced Course of the Swiss Society of Astrophysics and Astronomy, eds. A. Benz, 
M. Huber and M. Mayor, [Geneva Observatory: Sauverny], p. 3.). 
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Fellgett described the first numerically transformed two-beam interferogram and 
applied the multiplex method to stellar spectroscopy (Fellgett, P., J. de Physique et le 
Radium V. 19, 187, 236, 1958). Fellgett employed a Michelson-type interferometer 10 as 
shown in Figure 1, wherein the incoming beam of light B is divided into two beams Bi 
and B2.by a beamsplitter ("beam divider") 12. Bi is reflected from retro reflector 14, 
while B2 is reflected from retroreflector 16. As shown, beams Bi and B2 follow separate 
paths whose lengths can be precisely adjusted by delay lines (DLs) established by 
repositioning one or both of reflectors 14 and 16, as shown with reflector 16 connected to 
drive train 18 and drive motor 20. The beams, now with a path difference x (i.e. the 
"lags"), are recombined at beamsplitter 12 and focused by a concave mirror 22 onto a 
detector 24, producing an interferogram I(x), where 



as is discussed in more detail below. A mirror 28 is provided for an optional reference 
beam indicated in Figure 1 by the dotted lines. The reference beam is divided into two 
beams at the beamsplitter 12. These two beams are reflected from the retroreflectors 14 
and 16, are recombined at 12 and focused by 22 onto the detector 24. The reference beam 
allows the user to align the optics as well as to determine the zero-path position of the 
DLs, i.e., the position of retroreflectors 14 and 16 for which the path difference x = 0. In 
this manner, the reference beam thereby measures the delay, that is, the path difference x, 
introduced by the delay lines. This provides a more accurate determination of the optical 
path differences in the interferometer, and is typically included in applications involving 
the precise determination of spectral lineshapes and Doppler shifts. The detector is output 
to an amplifier and demodulator 30, and then the interferogram corresponding to the input 
spectrum is output to a recorder 32. 

The intensity of the combined beam is measured for a series of delay line positions. 
The wavelengths in the light beam cover a range from A, m j n to ^ maX5 i.e., centered on Xq and 

covering a range AX =A, max -A, m i n . The most important length parameter in the FTS is the 
lag, x, which is equal to the path length difference A-B. At any given wavelength X, 
complete constructive interference between light from the two paths occurs when x/X is an 
integer, and complete destructive interference occurs when x/X is an integer plus 1/2. 
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5 When the paths A and B are precisely equal to within a small fraction of Xq (i.e., x 

= 0 is the only delay for which x/X=0 at all wavelengths), the light waves at all 
wavelengths in the two beams constructively interfere and the intensity I in the 
recombined beam is at its maximum, I max . This position is known as the central fringe. As 
the DLs are moved and x changes, constructive interference between light waves from the 

1 0 two paths weakens, particularly at the shorter wavelengths, and I decreases. As the 
magnitude of x continues to increase, I reaches a minimum at x/Ao=l/2 and then rises 
again to a new (but weaker) maximum at x/Xq=1 . This weakening oscillation of I 
continues as x increases. When x/Xo is increased to many times Ao/AX, some wavelengths 
interfere constructively and some destructively, so I is close to the mean light level. Thus, 

1 5 if the observed spectral region is wide, there is only a small range of lag with large 
deviations from the mean level. 

The resulting data set of intensity measurements, I(x), measured at many values of 
x is known as an interferogram (Equation 1) as discussed above. The region of x over 
which there are large deviations from the mean level is termed the fringe packet. An 

20 example of a typical interferogram is shown in Figure 2. The wavelength of the high 

frequency oscillations is the central wavelength of the bandpass, X. As shown, the number 
of fringes in the central fringe packet is approximately equal to X/AX, where AX is the 
bandwidth. 

Typically, the interferogram is sampled in steps of Xo/2, and is then Fourier 
25 transformed to produce a spectrum. The spectrum is given as a series of values at regularly 
spaced discrete values of the wavelength, X. The spectrum that results from the Fourier 
transform of the interferogram contains artifacts of the PSF, which results from the finite 
lag range and the actual sampling of the interferometer. A wide range of deconvolution 
methods have been developed to disentangle the real signal from the deleterious effects of 
30 sampling, noise, etc., and have generally done so by implicitly modeling the spectrum as 
differing from zero only at discrete values of wavelength. The disadvantage of the 
deconvolution approaches is that they are highly nonlinear processes, so their behavior and 
uncertainties are hard to understand quantitatively. In addition, the disadvantage of 
modeling the spectrum only at discrete points is that the corresponding interferogram has 
35 significant sidelobes. 

The resolution of the spectrum at a given wavelength, X, is determined by the 
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5 maximum value of x/X and can be understood as follows. The light waves that comprise a 
narrow spectral line occupy only a small range of wavelengths, and thus stay correlated for 
a relatively long time, given roughly by bX/c, where 5X is the range of wavelengths 
making up the line and c is the speed of light. Since the lag x corresponds to a time delay 
between beams of x/c, a narrow line produces interference fringes over a large range of x. 

10 The FTS can measure over only a finite range in x, so it cannot distinguish between a 

spectral line of width SX and a narrower line that produces fringes over a larger range of x. 
For a spectral line wider than the resolution of the FTS, the width of a spectral feature is 
measured by the range of x over which there are interference fringes. 

The most common type of spectrometer is a dispersing spectrometer, consisting of 

1 5 a dispersing element (usually a grating) and a camera equipped with an array of detectors 
(usually a CCD) for multiplexing the dispersed output. Present CCD designs allow the 
number of channels N ch to exceed several thousand, so that the entire integration time is 
directed to integrating on all N ch channels. Recent planetary detections have used 
dispersing spectrometers with an absorption cell positioned in the path of the incoming 

20 beam to impose a reference set of spectral lines of known wavelength on the stellar 
spectrum. 

In principle, an FTS offers at least three major advantages over a dispersing 
spectrometer. First, the spectral resolution can be changed simply by changing the 
maximum value of the lag; second, the wavelength scale in the resulting spectrum is 

25 determined only by the delay line settings, while remaining insensitive to such effects as 
scattered light and flexure of the instrument; and third, the point spread function (PSF) of 
the spectrum can be determined to a high degree of precision. 

An FTS, however, also suffers certain disadvantages. These include low 
sensitivity: a conventional FTS is essentially a single-pixel scanning interferometer, and 

30 high spectral resolution requires measurements at a large number of lag settings. 

Accordingly, FTSs are commonly used when sensitivity is not a paramount concern, such 
as with laboratory spectroscopy or solar observations, or when very high spectral 
resolution or accurate wavelength calibration is required, such as in observations of 
planetary atmospheres. Other applications of FTSs include FTIR, MRI, and fluorescence 

35 and Raman emission spectroscopy. 

It would therefore be desirable to provide a spectrometer which offers the 
advantages of an FTS spectrometer while preserving most of the sensitivity of a dispersing 
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spectrometer. It would also be desirable to provide an improved algorithm for recovering 
the spectrum from the interferogram with greater fidelity, with easily quantifiable error 
estimates, and without producing undesirable artifacts. 



1 0 Summary of the Invention 

According to the invention, a dispersing Fourier Transform Spectrometer (DFTS) 
interferometer includes a Fourier Transform Spectrometer having an input for receiving a 
source light and an output, and a dispersive element having an input coupled to the Fourier 
1 5 Transform Spectrometer output and an output for providing the resulting multiple 

narrowband interferogram outputs of different wavelengths representative of the source 
light input. 

Also according to the invention, a method of determining a spectrum of a light 
source includes receiving and collimating a source light along a first optical path; 

20 transmitting a first part of the collimated source light further along the first optical path 
while reflecting a second part of the collimated source light along a second optical path; 
reflecting back the first part of the collimated source light along the first optical path; 
reflecting back the second part of the collimated source light along the second optical 
path; introducing a path length difference x between the first and second optical paths; 

25 recombining the back-reflected first and second parts of the collimated source light; 
dispersing the recombined beam into a plurality of different wavelengths; separately 
sensing an intensity I of each of the plurality of different wavelengths to thereby produce a 
set of data of interferogram intensities 1^ measured at a set of discrete lags Xj; and 
processing the data so as to produce a spectral output having a best fit with the set of data. 

30 The data processing preferably includes applying a sparse sampling algorithm for 

determining the best fit between a model interferogram and the data interferogram. In one 
form, the sparse sampling algorithm processes the set of data interferograms, Id(xi), where: 

d$ J t (s) C08(2TTXiS) 1 

35 s is the wavenumber, equal to the inverse of the wavelength, J t (s) is the true spectral 

intensity at wavenumber s, and the subscript / indicates that J t (s) is the truth spectrum and 
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5 is an unknown, and the wavenumbers Smi n (n) and s max (n) span the range of wavenumbers 
detected by the n* member of said set of light intensity sensing elements. A model 
spectrum, J m (sj), is selected, from which is inferred a model interferogram specified at a 
discrete set of lags x;, I m (xj); and a difference between the model interferogram and the 
data interferogram is determined, and an optimization method applied to determine a 

10 model interferogram best matched to the data interferogram Id(xj). 

The DFTS interferometer improves the sensitivity of a standard FTS by including a 
dispersive element, increasing the SNR by a factor of (Rg) 1/2 as compared to the FTS, 
where R g is the resolving power of the conventional dispersing spectrometer (i.e. R g = 
7JAX). The DFTS interferometer disperses the recombined light beam from the FTS 

1 5 module and focuses it onto a CCD detector, essentially splitting a single broadband FTS 
into N C h parallel, narrowband channels. A narrowband beam yields a spectrum with a 
higher signal-to-noise (SNR) ratio than a broadband beam because a narrowband beam 
filters out noise from wavelengths outside the bandpass without losing signal from inside 
the bandpass. The noise in the spectrum is a constant with a value proportional to the 

20 square root of the mean flux level in the entire interferogram, and restricting the bandpass 
and recording a narrowband interferogram filters noise from the spectrum without losing 
signal. 

The DFTS interferometer obtains spectra over a wide bandpass, with an easily 

configurable spectral resolution that can be very high (X/AX>\0^X with high sensitivity 
25 (i.e., high spectral SNR), with a well-known PSF, and with high velocity/wavelength 
precision (5A,M,=5v/c«10" 9 using a standard frequency stabilized HeNe laser for 
metrology). The algorithm is a preferred embodiment that models the spectrum as a 
continuous function rather than as a series of infinitely narrow delta functions as is done in 
conventional spectral deconvolution, and it solves the forward problem, i.e., it selects the 
30 set of spectral intensities J m (Sj) that yields an interferogram I m (xi) that best matches the 
measured interferogram Id(xi). 

The combination of the DFTS interferometer and the SSA provides more precise 
spectral measurements, a more precise position of the central fringe, and greater 
information about noise in the data. 
35 Astronomical spectroscopy is one application of the DFTS interferometer. The 

detection of planets requires the acquisition of sensitive, high-resolution, high-stability 
spectra from their parent stars. The minute oscillations in the Doppler velocity of the star 
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5 due to an orbiting planet have presented the most abundant signature of extrasolar planets 
to date. Jupiter-mass planets typically cause the lines in solar-type stars to be Doppler 
shifted by AA, sh j ft A«10" 8 . The detection of earth-mass planets is made possible with the 
DFTS interferometer when precisions improve to better than 3 m/s. 

The scope of other applications is broad. Spectroscopy is most commonly used to 
10 determine the composition of a sample. The desired spectral features are often faint and 
appear over a broad spectral range. The DFTS interferometer is ideal for such applications 
since the desired spectral regions can be isolated for maximum sensitivity. Precise 
knowledge of the PSF afforded by the DFTS interferometer can permit the identification 
of weak features juxtaposed to or superposed on intense background features, a capability 
1 5 owing to the FTS instrumentation component. 

Additional applications of the invention include but are not limited to manufacturing 
and product quality control, field detectors for hazardous compounds, Raman 
spectroscopy, radar detection algorithms, and atmospheric and metallurgic spectroscopy. 
Additional features and advantages of the present invention will be set forth in, or 
20 be apparent from, the detailed description of preferred embodiments which follows. 

Brief Description of the Drawings 

25 FIG. 1 is a schematic diagram of a prior art interferometer. 

FIG. 2 is a graph showing a representative conventional interferogram. 

FIG. 3 is a schematic diagram of a DFTS interferometer according to the invention. 

30 

FIG. 4 is a spectrum showing a representative FFT applied to a sparsely sampled 
interferogram. 

FIG. 5 is a spectrum showing the results of applying the SSA to the same dataset as in 
35 FIG. 4 according to the invention. 
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5 Detailed Description of the Invention 

Definition: The term "optics" as used herein when referring to a component of an 
interferometer of the invention includes a lens or a mirror. 

Referring now to Figure 3, a DFTS interferometer 100 includes a beamsplitter 102 

10 that partially reflects and partially transmits an input light beam B, splitting it into 
reflected beam Bi along a first optical path and transmitted beam B 2 along a second 
optical path. Bi is reflected from a first retroreflector 104 back to the beamsplitter 102, 
while B 2 is reflected from a second retroreflector 106 back to beamsplitter 102. As 
discussed above with respect to Figure 1, the paths of either or both of Bi and B 2 can be 

1 5 adjusted, such as is shown with respect to B 2 where a programmable drive-train such as 
that illustrated in Figure 1 is coupled to reflector 106 and thus introduce a path difference 
x between the first and second optical paths. Beams Bi and B 2 recombine at beamsplitter 
102 and due to path difference x produce an interferogram the properties of which are 
dependent on the spectral content of the incident input light beam B and on the optical 

20 path difference x. 

Unlike a conventional FTS that would then focus the light on a detector, the light 
from recombined beams Bi and B 2 is then directed onto a dispersing grating 108 to 
separate the beam into a plurality of channels differing in wavelength as is further 
described below. As illustrated in Figure 3 for a single channel, the light from each 

25 narrowband channel is focused by a lens 110 onto a detector, CCD 112. By dispersing the 
recombined light beam from the FTS module and focusing it onto a CCD detector, the 
DFTS interferometer splits a single broadband FTS output into N ch parallel, narrowband 
channels. A narrowband beam yields a spectrum with a higher signal-to-noise (SNR) ratio 
than a broadband beam because, unlike a conventional FTS spectrometer, the noise at any 

30 given wavelength in the spectrum is not proportional to the square root of the signal level. 
The noise in the spectrum is a constant with a value proportional to the square root of the 
mean flux level in the entire interferogram, and restricting the bandpass and recording a 
narrowband interferogram filters noise from the spectrum without losing signal. The 
DFTS interferometer realizes an increase in the SNR by a factor of (R g ) 1/2 as compared to 

35 a conventional FTS. This can be shown as follows. Consider a telescope collecting a 
stellar flux of W photons s' 1 nm" 1 (we have expressed W using units of wavenumber 
instead of wavelength). An interferogram with measurements at Ni ag delays is obtained 
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5 with a mean level of W ti ag As photons per lag for a given spectral channel, where As = 
s/Rg is the bandwidth of the channel, and t) ag is the integration time at each delay. In the 
following analysis, we consider the data from a single spectral channel Since the integral 
of the spectral intensities over the total spectral bandwidth is equal to the intensity, I Q , at 
the peak of the central fringe of the interferogram the mean spectral intensity (i.e., the 
10 mean signal level of the spectrum) is just I 0 divided by the spectral bandwidth. Assuming 
that the fringe contrast is 100%, then I 0 is just equal to the mean level of the interferogram, 
and the mean spectral intensity is: 

(3) 

On average, the noise level in the interferogram is determined according to Poisson 
1 5 statistics: 

(4) 

Parceval's Theorem states that the total noise power in the spectral and lag domains is 
equal: 



As 

(5) 

20 where a s is the average spectral noise power per pixel, and Oi is the average noise power 
in the interferogram per pixel. We combine the above equations to compute the signal-to 
noise ratio in the spectrum: 



SNR S = M^E. 

V #FTS m 

Not surprisingly, the number of samples in the interferogram (Ni ag ), is directly 
25 proportional to the number of independent spectral values, M, across one channel: 

M _ RftS _ 2iViag 

* f ' (7) 

and Equation 6 becomes: 
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SNR S = 




2Wti &g N\ iig s 



(8) 



For the case of the conventional FTS, R g = 1 . 

Since the width of the central fringe packet is inversely proportional to M, small 
values of M mean that meaningful signal is collected throughout a larger portion of the 
interferogram. A large value of M suggests a narrow fringe will be the only region in the 
10 interferogram that has significant signal. In effect, M serves to dilute the signal as the 
fringes decorrelate. Equation 8 demonstrates that SNRs is directly proportional to (R g ) 1/2 
for a constant integration time (t] ag Ni ag ), source brightness (W), observing wavenumber 
(s), and spectral resolving power (Rfts). Sensitivity is gained with greater multiplexing. 



1 5 includes an optional metrology detector for determining the path difference x. Two light 
beams enter the DFTS interferometer, one from the source to be measured ("science 
light") and the other from the laser metrology system ("metrology light"). Laser light from . 
the metrology laser is split into two beams with orthogonal polarizations at BSC1 . The two 
beams are frequency shifted (AOMs), recombined (BSC2), and spatially filtered and 

20 expanded (SPF) to the same size as the science light beam. Part of the recombined beam 
is split from the main beam (B), both orthogonal polarizations are mixed at the polarizer 
(PI) and focused onto a reference detector (Dl). Light from the source to be measured 
enters the spectrometer through a polarizer (P3). Polarized science light is combined with 
metrology light at the notch filter (Nl). The combined beam is split into two by a 

25 polarizing beamsplitter (BSC3). Each beam propagates through DL1 or DL2. The beams 
are recombined at a polarizing beamsplitter (BSC4). The metrology light is separated 
from the combined light using a notch filter (N2), the orthogonal polarizations are mixed 
with a polarizer (P2) and sent to the metrology detector (D2). The intensity measured at 
D2 is compared with that measured at Dl to generate the metrology signal. At this point, 

30 a conventional FTS would focus the light transmitted through N2 to a detector. Instead, 
with the DFTS interferometer, the light is sent to a dispersing spectrometer as shown. The 
two polarizations are separated with a Wollaston prism (W), dispersed with a transmission 
grating (G), and are each focused onto a row of pixels on the CCD 112.The data in the 
form of interferogram intensities, I d , measured at a set of discrete lags, x i3 where l<i<N, is 

35 recorded on the CCD 112 with a computer. 



Figure 3 illustrates a DFTS interferometer 200 according to the invention, that 
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The DFTS interferometer in a preferred embodiment utilizes the concept that the 
interferogram obtained may be "undersampled" as compared to the sampling required 
using a conventional FTS. According to the Nyquist Theorem, a FTS must be sampled at 
increments 5x of the path difference A-B such that 5x=l/(2As), where As=s max -s m j n is the 

width in wavenumbers of the spectral region being observed. (For a given wavelength X, 
the wavenumber s is l/X.) 

As an example, with a standard FTS the wavenumber range As st< j is large, so the 
sampling interval 5x stc j is small. The DFTS interferometer, however, functions as R g 
standard FTSs working in parallel (where R g is the number of spectral channels in the 
dispersing spectrometer), each working in a narrow band. The sampling interval 8x g for 
these narrow-band FTSs is l/(2As g ), where As g is the wavenumber range of a single 
channel. Because the channels are 1/R g as wide as the bandpass of the standard FTS, 
As g =As st( j/R g , and the sampling interval 5x g is R g times larger than 5x std . Therefore, the 

number of sampled points needed to attain a given resolution is reduced by a factor of up 
to R g for each channel. Also, conventional FTS data processing techniques involve 

converting an interferogram into a spectrum using a Fourier transform, after which 
corrections are applied for sampling to result in a final spectrum. The DFTS technique 
instead applies an algorithm for inferring the best spectrum given a bandwidth limited 
interferogram. 

The algorithm is preferably applied as either (1) a fast algorithm for solving the 
best spectrum assuming that the location of the central fringe for each channel is known, 
or (2) a slow algorithm, for solving the best spectrum as well as the best location for the 
central fringe. In both approaches, the data consist of interferogram intensities, 1^, 

measured at a set of discrete lags, x i? where l<i<N. Apart from a constant, which can be 

ignored in this analysis, the interferogram is simply the inverse cosine transform of the 
spectrum. Therefore, the data can be written as: 




(9) 
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5 where J t (s) is the spectral intensity at wavenumber s. The subscript t indicates that J t (s) is 
the truth spectrum, and is not known to the observers. It is then desired to infer J t (s) based 
on observations of I<j(Xj). 

At this point, there are two significant departures from conventional approaches. 
The first is that the forward problem is solved. The forward problem is the process of 

10 selecting the set of spectral intensities, J m (sj), which yields an interferogram I m (Xj) that 
best matches Id(xj). This is the opposite of the standard strategy of solving the backwards 
problem by doing a deconvolution of I^Xj) in the hopes of disentangling the real signal 
from the deleterious effect of sampling, noise, etc., and recovering J t (sj). 

The second departure from conventional methodologies is that a model spectrum 

1 5 having continuous frequency coverage is selected. Conventional methods apply Fourier 
Transforms to discretely sampled data and return discrete data. The results from 
conventional methods are diminished in quality due to the lack of knowledge between 
sampled frequencies. The method presented here alleviates this problem to first order. 

Initially, one starts by guessing a set of M spectral intensities, J m (Sj), which span a 

20 wavenumber range defined by the edge wavenumbers of a single, narrowband spectral 

channel. It is known that the light outside this wavenumber range has been excluded from 
the detector by the conventional spectrometer in the FTS optical train. Furthermore, it can 
be assumed that the continuous spectral intensities between sj and Sj+j are given by the 
interpolation between J m (sj) and J m (Sj +1 ). 

25 This choice of J m (sj) and the assumptions above result in an interferogram given 

by: 

M-l 

Im(Xi) = 22 I d 3 [Jm(Sj) + {s - Sj)Aj] COS (2^$), 

31 (10) 

where: 



30 



S j+l - s j 



(11) 

and 8 is the location of the central fringe in the interferogram. The integral can be 
evaluated analytically, reducing the expression to: 

12 
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M-l 



Im{Xi) = ^2 [<Xi,jJ m (8j) + 
j-1 



where: 



(12) 



CXij = 



sin (27rXjSj + i ) - sin (27r^Sj) 



and 



(13) 



ft,; = 



{sj+i — 8j)sm(2wXiSj+i) 



+ 



10 



' cos (2wXiSj+i) - cos(27ra?iSj) 

(14) 



In the fast algorithm technique, the variance of the residuals between the model 
interferogram and the data interferogram is given by: 

n 



(15) 



1 5 Above, when describing the forward problem, it is desired to obtain a model 

interferogram best matched to the data interferogram. This condition can be expressed as 
a set of equations: 



dr 



dJ m (s 



To complete the problem requires the Jacobian, which can be derived analytically. 



(16) 



20 



\dJm(Sj)J V 5 2-Sl/ 



i = i, 



(17) 



= (_&!=!_) + ^ _ for 2 < > < M - 1, 



(18) 



25 



\dJ m {Sj)J \sm - sm-i J 



(19) 
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The slow algorithm technique also starts with the expression for Equation 15. 
Equation 16 is still valid, and in addition: 

!Sa-!±w*-«- w («*4)-. 

<=1 (20) 
The Jacobian is given by Equations 17 — 19 and: 

> =1 (21) 

where: 

sinfa^+i) smjzjSj) 



A i,3 = ~ S j+l COS(Zi$ j+1 ) + Sj COS{Zi$j) + 



2i Zi 



15 (22) 
and: 

= Sjgj+lCOSfes i+ i)+(23^^ 

cos(, iS , +l ) + ^^-^^, 

Z i Z i (23) 

20 We have used the definition zj=2n(xj-z). 

Accordingly, with the spectral reconstruction algorithm described above, a model 
spectrum is first selected, resulting in a model interferogram. The model spectrum is then 
varied to yield a model interferogram that most closely matches the data interferogram in a 
least-squares sense, a method that may be implemented using a simple Newton technique. 

25 The algorithm functions primarily as an anti-aliasing filter, replacing the traditional 

Fast Fourier transform (FFT) for our application. The spectrum resulting from a sparsely 
sampled interferogram consists of the "true" spectrum plus an "aliased" version of the 
spectrum shifted to other wavelengths. So long as the sampling interval in the 
interferogram satisfies the Nyquist Theorem, the aliases will not overlap the "true" 

30 spectrum. The algorithm infers the spectral intensities only over a user-defined bandwidth 
of interest (presumably containing the "real" signal and not one of the aliases). These 
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effects are illustrated in Figure 4, which shows the results of applying a FFT to a sparsely 
sampled interferogram. The real spectrum (denoted by the arrow) is aliased at all 
frequencies. As the sampling in the interferogram approaches the Nyquist Limit, that is 
becomes sparser, the aliased peaks merge. Figure 5 shows the results of applying the 
algorithm according to the invention for the same dataset as used with the FFT in 
generating Figure 5. It is evident from comparing the figures that the sparse 
reconstruction algorithm serves as an anti-aliasing filter, and more efficiently reconstructs 
the actual signal than does the FFT. The spectrum outside this interval is zero. 

The DFTS interferometer of the invention is an achromatic device in that it can 
obtain spectra from a luminous source at any given wavelengths throughout the 
electromagnetic spectrum, so long as the optical components and detectors are selected so 
as to provide reasonable sensitivity. Obviously many other modifications and variations 
of the present invention are possible in the light of the above teachings. It is therefore to 
be understood that the scope of the invention should be determined by referring to the 
following appended claims. 
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